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Abstract 

Discrete breathers are time-periodic, spatially localized solutions of equa¬ 
tions of motion for classical degrees of freedom interacting on a lattice. They 
come in one-parameter families. We report on studies of energy properties 
of breather families in one-, two- and three-dimensional lattices. We show 
that breather energies have a positive lower bound if the lattice dimension of 
a given nonlinear lattice is greater than or equal to a certain critical value. 
These findings could be important for the experimental detection of discrete 
breathers. 
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Recently progress has been achieved in the understanding of localized excitations in 
nonlinear lattices. Discrete breathers (DBs) are time-periodic, spatially localized solutions 
of equations of motion for classical degrees of freedom interacting on a lattice [|1|], [Q, 
Nowadays it is known that the reason for the generic existence of DBs is the discreteness of 
the system paired with the nonlinearity of the differential equations defining the evolution 
of the system [£f|, ||. Thus one can avoid resonances of multiples of the discrete breather’s 
frequency fib with the phonon spectrum fl q of the system If the coupling is weak the 
phonon spectrum consists of narrow bands. The nonlinearity and the narrowness of the 
phonon bands allows for periodic orbits whose frequency and all its harmonics lie outside 
the phonon spectrum. For some classes of system, existence proofs of breather solutions 
have been published [0, ||, ||. A list of references is given in 

For generic Hamiltonian systems, periodic orbits occur in one-parameter families, and 
discrete breathers are no exception. In many cases, the energy can be used as parameter 
along the family, but as is well known, the energy can have turning points along a family 
of periodic orbits. Mathematically, such a turning point in energy is called a saddle-centre 
periodic orbit. 

The main message of this paper is that in 3D lattices, a turning point (in fact, minimum) 
in energy is almost inevitable for discrete breather families. 

One important property of DBs is their generic existence for weak enough coupling, 
independent of the lattice dimension [?], |7|]. This means that DBs are not just a ID curiosity 
but could be interesting from the point of view of applications. The experimental detection 
of DBs requires some additional knowledge about their properties. In this contribution we 
give heuristic arguments that the energy of a DB family has a positive lower bound for 
lattice dimension d greater than or equal to some d c , whereas for d < d c the energy goes 
to zero as the amplitude goes to zero, and we confirm these predictions numerically. The 
critical dimension d c depends on details of the system but is typically 2 and never greater 
than 2. Furthermore, for d > d c , the minimum in energy occurs at positive amplitude and 
finite localisation length. Consequently experiments could be designed to look for activation 
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energy thresholds for localized excitations. 

Let us consider a d-dimensional hypercubic lattice with N sites. Each site is labeled by 
a d-dimensional vector l e Z d . Assign to each lattice site a state Xi e Rf, where / is the 
number of components and is to be finite. The evolution of the system is assumed to be 
given by a Hamiltonian of the form 

H = J2 H loc {X{) + H int {Xu £W), (1) 

i 

where H int depends on the state at site l and the states X/ /+s in a neighbourhood. We assume 
that H has an equilibrium point at Xi = 0, with H({Xi = 0}) = 0. 

DB solutions come in one-parameter families. The parameter can be the amplitude 
(measured at the site with maximum amplitude), the energy E or the breather frequency 
fib- It is anticipated (and was found both numerically and through some reasonable ap¬ 
proximations [|]) that the amplitude can be lowered to arbitrarily small values, at least for 
some of the families for an infinite lattice. In this zero amplitude limit, the DB frequency 
f \ approaches an edge of the phonon spectrum fl q . This happens because the nonresonance 
condition fl q /flb ^ 0,1, 2, 3,... has to hold for all solutions of a generic DB family p]. In the 
limit of zero amplitude, the solutions have to approach solutions of the linearized equations 
of motion, thus the frequency fib has to approach some fl q , but at the same time not to 
coincide with any phonon frequency. This is possible only if the breather’s frequency tends 
to an edge f 1^ of the phonon spectrum in the limit of zero breather amplitude. If we consider 
the family of nonlinear plane waves which yields the corresponding band edge plane wave 
in the limit of zero amplitude A, then its frequency fl will depend on A like 


|D — flj 


A z 


( 2 ) 


for small A, where the “detuning exponent” z depends on the type of nonlinearity of the 
Hamiltonian ([[]), and can be calculated using standard perturbation theory HIT] . 

It is tempting to check then whether the breather appears through a bifurcation from 
a periodic orbit which is a normal mode of the linearized equations of motion for any 
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system with finite N. Band edge plane waves of the linearized equations of motion can be 
continued to non-zero amplitudes for the general nonlinear system. The stability analysis 
of these periodic orbits yields the possibility of tangent bifurcations (collision of Floquet 


multipliers at +1) if some algebraic inequalities of the expansion coefficients of H in (|I) are 
met . It has been also shown that the orbits which bifurcate from the plane wave are not 
invariant under discrete translations and have the shape of discrete breathers [jl2[]. It has 
been conjectured that the new bifurcating orbits are discrete breathers. Subsequently it was 
successfully explained why discrete breathers exist or not for certain models by analyzing 
the above-mentioned algebraic inequalities [12[]. Numerical studies confirm these findings 
13| for some one-dimensional models. 

The above-mentioned analysis of stability of band edge plane waves was carried out for 
systems with detuning exponent z — 2 and large N. The critical amplitude A c of the plane 
waves at the bifurcation point depends on the number of lattice sites as A c ~ N~ l ^ d |L2| . 


We see that the amplitudes of the new orbits bifurcating from the plane wave become small 
in the limit of large system size. If the energy of the system is given by a positive definite 
quadratic form in the variables X in the limit of small values of X it follows for the critical 
energy of the plane wave at the bifurcation point W% 


E c ~ N 


1-2 /d 


( 3 ) 


Result (|3|) is surprising, since it predicts that for z = 2 the energy of a DB for small 
amplitudes should diverge for an infinite lattice with d — 3 and stay finite (nonzero) for 
d = 2, whereas if d — 1 the breather energy will tend to zero (as initially expected) in 
the limit of small amplitudes and large system size. The whole construction depends on 
the validity of the assumption that the new periodic orbits bifurcating from the plane wave 
through the above-mentioned tangent bifurcation are indeed DBs. 

It is not known how to prove this assumption. But we can estimate the discrete breather 
energy in the limit of small amplitudes and compare the result with (3j). Define the amplitude 
of a DB to be the largest of the amplitudes of the oscillations over the lattice. Denote it by 
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Aq where we define the site / = 0 to be the one with the largest amplitude. The amplitudes 
decay in space away from the breather center, and by linearising about the equilibrium state 
and making a continuum approximation, the decay is found to be given by Ai ~ CFd(\l\S) 
for |/| large, where F ( j is a dimension-dependent function 


Fi(x) 



d( , 


( 4 ) 

( 5 ) 


5 is a spatial decay exponent to be discussed shortly, and C is a constant which we shall 
assume can be taken of order A 0 . To estimate the dependence of the spatial decay exponent 
5 on the frequency of the time-periodic motion fife (which is close to the edge of the linear 
spectrum) it is enough to consider the dependence of the frequency of the phonon spectrum 
fl q on the wave vector q when close to the edge. Generically this dependence is quadratic 
(flE — fig) ~ |q — qE \ 2 where fi^ ^ 0 marks the frequency of the edge of the linear spectrum 
and qE is the corresponding edge wave vector. Then analytical continuation of (q — g#) to 
i (q — qE) yields a quadratic dependence [fife — fi^| ~ S 2 . Finally we must insert the way 
that the detuning of the breather frequency from the edge of the linear spectrum [fife — fi#| 
depends on the small breather amplitude. Assuming that the the weakly localized breather 
frequency detunes with amplitude as the weakly nonlinear band edge plane wave frequency 
this is | fit — fi#| ~ Aq. Then S ~ Aq 2 . 

Now we are able to calculate the scaling of the energy of the discrete breather as its 
amplitude goes to zero by replacing the sum over the lattice sites by an integral 


E b ~ 



r d 1 Fj(5r)dr ~ 


^(4-*d)/2 


( 6 ) 


This is possible if the breather persists for small amplitudes and is slowly varying in space. 
We find that \i d > d c = A/z the breather energy diverges for small amplitudes, whereas 
for d < d c the DB energy tends to zero with the amplitude. Inserting z — 2 we obtain 
d c = 2, which is in accord with the exact results on the plane wave stability [12] and thus 
strengthens the conjecture that discrete breathers bifurcate through tangent bifurcations 
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from band edge plane waves. Note that for d = d c logarithmic corrections may apply to ([|), 
which can lead to additional variations of the energy for small amplitudes. 

An immediate consequence is that if d > d c , the energy of a breather is bounded away 
from zero. This is because for any non-zero amplitude the breather energy can not be zero, 
and as the amplitude goes to zero the energy goes to a positive limit (d = d c ) or diverges 
(d > d c ). Thus we obtain an energy threshold for the creation of DBs for d > d c . This 
new energy scale is set by combinations of the expansion coefficients in (|l ). If z = 2 with 
|0 — Qe| ~ (3A 2 for the nonlinear plane waves, and the energy per oscillator E ~ gA 2 and 
the spatial decay exponent 6 is related by ~ k5 2 1 then the energy threshold E min 

is of the order of K,g//3, and the minimum energy breather in 3D has spatial size of the 
order of the lattice spacing, independently of k, g and f3. One should allow for a factor of 
(2 + d) for underestimating the true height of the minimum and the contributions of nearest 
neighbours. 

To confirm our findings, we performed numerical calculations. First we study the discrete 
nonlinear Schrodinger (DNLS) equation 

T^i^ + IV^ + C £ *m), (7) 

m^Ni 

where Ni denotes the set of nearest neighbours of l. The detuning exponent z is easily 
seen to be /r — 1. Making the substitution dq = Aie lQbt we solve the algebraic equations for 
the real amplitudes Ai. Numerically this is implemented by considering the case of large 
breather amplitude A 0 first. Then the breather is essentially given by A 0 ps (fit — l) 1 / 
and Ai^o = 0. Next we define a functional G which is the sum over the squares of differences 
between left hand and right hand parts of all algebraic equations for the amplitudes. This 
functional is minimized by gradient descent, where the initial guess is the large amplitude 
approximate solution. Finally the frequency flj, is varied in small steps and the breather 
solution is traced. In Fig.l we show the resulting breather energy as a function of the 
amplitude Aq for /j, = 3 and d = 1,2,3. The results are in full accord with the predictions. 
For d = 3 the above estimate of the minimum energy yields a value of 0.2 with (3 — g — 1 
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and k — C — 0.1. The mentioned factor (2 + d) = 5 accounts for the deviation from 
the true value of 1. Fig.2 shows the amplitude distribution of the discrete breather with 
minimum energy in the (x, y)-plane crossing the breather center for d = 3. The minimum 
energy breather is strongly localized - its spatial width is only few lattice spacings. In Fig.3 
we show results for d — 1 and fi — 3, 5, 7. Again we find full agreement. Note that even 
one-dimensional lattices exhibit positive lower bounds on breather energies if n > 5. This 
d = 1 result has also been predicted using variational techniques |L4 . 


To demonstrate that the numerical results are not an artefact of the DNLS case, we 
study the d -dimensional nonlinear Klcin-Gordon lattice 


Ui = -Ui -U?-C £ (Ut - U m ). ( 8 ) 

mGNi 

The detuning exponent z is given by /i — l for /i odd and 2/i — 2 for /i even. Again the discrete 
breather with large amplitude is essentially an on-site excitation and given by Uq = —Uq — Uq 
and Uty o = 0. The equations of motion are integrated numerically for a given set of initial 
conditions {Ui(t = 0),Ui(t = 0)} over the breather period T b = 27r/Q b . The functional G = 

(( Ui(T b ) — Ui( 0)) 2 + ( Ui(T b ) — f/;(0)) 2 ^ is minimized with respect to the initial conditions 
using gradient descent. This method allows us to perform a reliable numerical calculation 
of DBs in 3-dimensional arbitrary lattices. The result in Fig.4 for fi = 3 and d — 3 is again 
in full accord with the predictions. 

We can predict that a modified DNLS system with an additional term ty |4q| Al, ~ 1 'Ih can 
exhibit complex curves E b (A 0 ). For example, for d — 1, fi — 7, fi' — 3 and ty = 0.1, the 
E b [A {) )-dependence will be nearly identical to the case ry = 0 already considered, if the 
amplitude Aq is not too small. Then E b (Ao) will show a minimum at a non-zero value of 
A 0 . For small A 0 however the energy of the breather will ultimately decay to zero, so the 
curve has a maximum for smaller amplitudes! The dashed line in Fig.3 shows the numerical 
calculation, which coincides with our prediction. 

Our findings should help to detect discrete breathers in experimental realizations like 
the dynamics of atoms in crystals. For a 3-dimensional crystal we predict a positive energy 
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threshold for the excitation of discrete breathers. 


Another consequence of our work is that breather solutions belonging to parts of the 
family where the energy is decreasing with increasing amplitude are dynamically unstable, 
whereas those in the other parts have a good chance of being dynamically stable. This can be 
seen from a Poincare map of the phase space flow around the breather orbits. The minimum 
energy breathers correspond to saddle-centre bifurcations, since no breather solutions will 
exist if the energy is lowered beyond the minimum breather energy. 

A similar phenomenon occurs in polaron theory. In a three-dimensional lattice, two 
polarons of unit electric charges exist above a certain parameter threshold (large and small 
polaron) ||T5|| . 

Summarizing, we have shown that discrete breather families have positive lower energy 
bounds if the dimension of the lattice is larger or equal to some critical value which in turn 
is defined by the power of the first nonlinear expansion term in the equations of motion. 
These results are expected to be of importance for the experimental detection of discrete 
breathers, because the minimum energy of a breather family should show up as an activation 
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FIGURE CAPTIONS 


Fig-1: 

Breather energy versus amplitude for the DNLS system in one, two and three lattice dimen¬ 
sions. Parameters C = 0.1 and n — 3 for all cases. System sizes for d = 1,2,3: 1V=100, 
1V=25 2 , A^=31 3 , respectively. The estimated points (A; E) of bifurcation of the band edge 
plane wave for d— 1,2,3 are: (0.014;0.024), (0.064;5.53), (0.097; 237), respectively. 

Fig-2: 

Amplitude distribution of the minimum energy breather solution of the DNLS system with 
d — 3, y — 3, C — 0.1 and N = 31 s . Actually only a distribution in a cutting (x;y) plane 
is shown (the plane cuts the center of the breather). The intersections of the grid lines 
correspond to the actual amplitudes, the rest of the grid lines are guides to the eye. 

Fig-3: 

Breather energy versus maximum amplitude for the DNLS system in one lattice dimension 
and for three different exponents fi = 3,5,7 (solid lines). The system size is IV = 100 and 
the parameter C — 0.1. The dashed line is for the modified system (cf. text). 

Fig.4: 

Breather energy E versus frequency detuning (D — D^) for a 3D Klcin-Gordon lattice. 
Parameters /i = 3 and C = 0.1. System size 1V=10 3 . 
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